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Disordered magnetism in 
superconducting KFe 2 As 2 single 
crystals 
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High-quality KFe2As2 (K122) single crystals synthe- 
sized by different techniques have been studied by mag- 
netization and specific heat (SH) measurements. The 
adopted phenomenological analysis of the normal state 
properties shows that there are two types of samples both 
affected by disordered magnetic phases: (i) cluster-glass 
(CG) like or (ii) Griffiths phase (G) like. For (i) at low 
applied magnetic fields the T-dependence of the zero 
field cooled (ZFC) linear susceptibility (xi) exhibits an 
anomaly with an irreversible behavior in ZFC and field 
cooled (FC) data. This anomaly is related to the freez- 
ing temperature Tf of a CG. For the investigated sam- 
ples the extrapolated Tf to B — varies between 50 K 
and 90 K. Below Tf we observed a magnetic hysteresis 
in the field dependence of the isothermal magnetization 
(M(B)). The frequency shift of the freezing temperature 
ST f = ATf/[T f A(\TLu)\ ~ 0.05 has an intermediate 
value, which provides evidence for the formation of a 



CG-like state in the K122 samples of type (i). The 
frequency dependence of their Tf follows a conven- 
tional power-law divergence of critical slowing down 

r = T [Tf(y)/Tf(Q) - l]- zv> with the critical expo- 
nent zv =10(2) and a relatively long characteristic time 
constant to=6.9- 10~ n s also supporting a CG behavior. 
The large value of the Sommerfeld coefficient obtained 
from SH measurements of these samples was related to 
magnetic contribution from a CG. Samples from (ii) did 
not show a hysteresis behavior for xi(T) and M (B). Be- 
low some crossover temperature T* ~40K a power-law 
dependence in the \i ^ T G , with a non-universal 
A G was observed, suggesting a quantum G-like behavior. 
In this case xi an d M (B) can be scaled using the scal- 
ing function M S (T,B) = B 1 - Xa Y( f iB/k b T) with the 
scaling moment /i ~3.5/x&. The same non-universal ex- 
ponent was found also in SH measurements, where the 
magnetic contribution C/T oc T Xg . 



X 



1 Introduction The interplay between superconduc- 
tivity (SC) and magnetism, and the role of correlation ef- 
fects in Fe-pnictides are under debate UEE]. ^ * s com- 
monly assumed that in most of the so-called stoichiomet- 
ric parent compounds nesting between electron (el) and 
hole (h) Fermi surfaces is responsible for the presence of 
long-range spin density waves (SDW). To get SC, the SDW 
should be suppressed by chemical doping or external pres- 
sure EEI- Therefore, it is believed that SC is driven at least 
partially by AFM spin fluctuations (sf). In contrast, in some 
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other of Fe-pnictides such as LiFeAs (111) or KFe2As2 
(K122) there is no nesting. Therefore, the role of mag- 
netism in the SC pairing of these compounds is less ob- 
vious but it is still believed that remaining sf as in Bal22 
||4l or a new type sf can be responsible for SC. 

In this paper we demonstrate that the physical prop- 
erties of clean K122 single crystals are strongly affected 
by an unexpected glassy like magnetic behavior such as 
spin-glass (SG) and Griffiths (G) type. It is known that a 
SG phase gives a nearly linear magnetic contribution to 
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Figure 1 (Color online) The temperature dependence of 
the volume susceptibility of SI and S2 in the SC state. In- 
set: Their molar susceptibilities in the normal state. 



the SH below the freezing temperature Tf [6 , 7 1 . In some 
cases the SG contribution can be hardly distinguished from 
the usual electronic (Sommerfeld) contribution to the SH, 
since it has the same linear T-dependence and only a weak 
maximum slightly above Tf. Therefore, the Sommerfeld 
coefficient and the deduced strength of correlation effects 
can be signifianfly overestimated, if one considers only SH 
data ignoring thereby the glassy phases. Moreover, the in- 
terplay of superconductivity (SC) and unknown magnetic 
phases can lead to confusing conclusions concerning SC 
gaps (H. A clear understanding of various coexisting or 
competing forms of magnetism should be addressed first. 

2 Experimental results and discussion 

2.1 Samples K122 single crystals have been grown 
using a self-flux method (FeAs-flux (SI) and KAs-flux 
(S2)). The high quality of the grown single crystals was 
assessed by complementary techniques. Several sam- 
ples were examined with a Scanning Electron Micro- 
scope (SEM Philips XL 30) equipped with an electron 
microprobe analyzer for a semi-quantitative elemen- 
tal analysis using the energy dispersive x-ray (EDX) 
mode (for details see (9 fT0l ). Resistivity of all mea- 
sured samples showing a metallic behavior at all T with 
a RRR 5 = p(300K)/p(5K) «400-500, where p(300K) 
and p(5 K) is resistivity at T = 300 and 5 K, respectively, 
comparable with the best values reported in the literature. 
The resistivity data will be published elsewhere. Low-T 
SH and ac susceptibility were determined using a PPMS 
(Quantum Design). The dc magnetic susceptibility has 
been measured in a SQUID (Quantum Design) magne- 
tometer. In this paper the data for two representative single 
crystals are shown. 

2.1.1 Magnetization measurements 

Samples with cluster glass behavior Fig. Q] de- 
picts the T-dependence of the volume susceptibility (x v ) 
determined from dc magnetization of our samples mea- 
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Figure 2 (Color online) (a) The molar susceptibility 
Xi(T) = M/B for SI crystal measured at different mag- 
netic fields B. M is the magnetization. Inset: field depen- 
dence of freezing temperature for field -Bii & and Bn c . (b) 
Field dependence of magnetization of S 1 crystal measured 
after ZFC at T=5K and T=60K (B|| o6 ) . Inset: field depen- 
dence of magnetization at Bu^ and Bu c and T=5K. 



sured under both zero-field-cooled (ZFC) and field-cooled 
(FC) conditions with the field B^ ab = 20 Oe. Bulk SC 
of our samples is confirmed by 'full' diamagnetic signals 
of the ZFC data at low T. For sample SI, a clear split- 
ting between ZFC and FC normal state linear susceptibil- 
ity Xi{T) = M{T)/B curves is observed below 100 K 
(see the inset of Fig. [TJ, where M (T) is a magnetization 
in the field B. The maximum in the ZFC Xi{T) is at- 
tributed to the freezing temperature of a spin glass (SG) 
type phase at T f « 65 K and B = 20 Oe 03] ■ T f de- 
creases with increasing field and at 5 T no splitting was 
observed down to 2 K (see Fig.|2^). The field dependence 
of Tf is shown in the inset of Fig. [2^. The extrapolated 
value of Tf(B = 0) ~ 90 K is a relatively high value. 
This might point to a large concentration of the involved 
magnetic moments (MM) >10% in sample SI Q. Such a 
high value of MM is expected from entropy estimations, 
too (see section 2.2). On the other hand, structural investi- 
gations did not reveal any impurity phase (see section 2.3). 
Therefore, we speculate that the high value of Tf might be 
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caused by a low-lying excited incommensurate spin den- 
sity wave state lfl2ll . For a more detailed consideration of 
this scenario see Ref. 13. In addition, an upshift of the max- 
imum and its lowering with increasing frequency v of the 
ac susceptibility, generic for a SG ordering (6][7) (Fig. |3), 
was observed for crystal SI. The value of the frequency 
shift of Tf HD: 



ST 



AT, 



f 



[T f A{logu)\ 



0.05. 



(1) 



is above the range 0.001-0.02 expected for canonical SG 
but well below ~ 0.3 observed in the case of superpara- 
magnets Q. Such an intermediate value of the frequency 
shift is usually related to a the so-called cluster glass (CG) 
behavior M14U15L The frequency dependence of the Tf 
shown in inset Fig.[3]follows a conventional power-law di- 
vergence of critical slowing down IT7l[T5l: 



To 



Tf(0) 



(2) 



where r = 1/v is the relaxation time corresponding to the 
measured frequency v, tq is the characteristic relaxation 
time of single spin flip, T f [y = 0,B = 50e) w71 K 
is the spin-glass temperature as the frequency tends to 
zero adopted from dc susceptibility measurements (inset 
Fig. |2k), and zv is the dynamic critical exponent. It is con- 
venient to rewrite Eq.[2]in the form: 



ln(r) = in(ro) — zv ln(t), 



(3) 



where t = Tj(v)/Tj(Q) - 1. The fit of the experimental 
data by a power-law divergence Eq.[3]is shown in the inset 
Fig. [3] The best fit was obtained with zv =10(2) and tq = 
6.9- 10~ u s. The value of zv is in the range of 4 to 12 
observed for typical SG fl31 . On the other hand the value 
of To is large as compared to 10 -12 — 10 _14 s observed for 
structureless SG systems, which is of the order of the spin- 
flip time of atomic magnetic moments (10 _13 s) fl6l . This 
suggests a slow spin dynamics in our crystal SI, likely due 
to the presence of interacting clusters rather than individual 
spins. 

Another signature of a SG-like behavior for crystal S 1 
is a hysteresis in the magnetization data below Tf with a 
reduced ZFC susceptibility at low fields (Fig.[2j)). This be- 
havior is expected in the case of SG or CG systems Q 
n~4lrT7Tl below Tf and also excludes superparamagnetic be- 
havior in our samples where no hysteresis was observed 
ifTHl . On the other hand, in the case of canted antiferromag- 
netic (AF) or ferromagnetic (FM) impurity phases hystere- 
sis is expected but with a higher susceptibility at low fields 
because the clusters are at first saturated along their local 
easy axis, and only after that various clusters become fully 
aligned along the applied field lfl6l . Therefore, our M(B) 
data exclude (large) clusters of impurity phases such as 
FeAs, Fe2As or other iron compounds. The same conclu- 
sion can be drawn from our SH measurements (see below). 
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Figure 3 (Color online) T-dependence of the real part 
of the ac susceptibility for sample S 1 measured for three 
different frequencies v at 5 Oe ac field amplitude. Inset: 
the zy-dependence of the freezing temperature plotted as 
ln{r) vs. ln(£), where t denotes the reduced temperature 
t = T f (v)/Tf(0) - 1. 



We also observed the displacement of the ZFC magnetiza- 
tion compared to magnetization after complete hysteresis 
loop at magnetic field applied parallel to ab plane (inset 
Fig-lit). For the field along the c-axis no displacement was 
observed. This indicates that the glassy behavior is domi- 
nated by a magnetic interaction between moments lying in 
the afe-plane. 

Samples with Griffiths-phase behavior In con- 
trast, the T-dependence of the linear susceptibility xi (^) 
of some crystals S2 does not show a difference between 
ZFC and FC curves above T c . The xi (T) data of one of the 
S2 crystals is shown in Fig. St- At high T >200 K, xi(T) 
follows a Curie- Weiss behavior with an AFM @ c =-l 17 K 
ED. At T <120 K xz(T) shows a plateau-like feature 
with a field independent susceptibility above B = 1 T. In 
a view of the observed CG in sample SI, we relate this 
flattening to a tendency to form magnetic clusters in the 
crystals S2, too. However, with lowering T after a weak 
reduction of xi {T) instead to form a CG phase in crystal 
S2, below T* an exponential increase of the susceptibility 
is observed: 



Xi(T) = Xia + 



Ci 



T 



1-A G 



(4) 



where xio is a T-independent susceptibility, and Cu a is 
a back ground constant. A power-law with the exponent 
A G was found up to the highest measured field 7 T (see 
Fig. [4^). This exponential behavior is quite similar to the 
reported one for the weak itinerant ferromagnetic alloy 
Nii-^Va; l20l . where the formation of a Griffiths (G) phase 
with a non-universal exponent was observed near a FM 
quantum critical point. Following the analysis proposed 
there, the field and T-dependence of the magnetization can 
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Figure 4 (Color online) (a) The molar susceptibility 
Xi(T) = MJB for crystal S2 measured in different 
magnetic fields B, where M is the magnetization. Fit- 
ting curves using Eq. (b) The scaled magnetization 
M S B~ - 5 = (M(T,B)~xwB)B-°- 5 vs. B/T for crystal 
S2.(seeEq.|5) 



be scaled on a single curve Fig . |4j^ : 

Ms (t,b) = b 1 -^yA 

to 2 



(5) 



where /i is the scaling moment and Y = A /(l + z 2 ) A °/ 2 
a scaling function with A = A/ /i Ac as a constant. To scale 
the data using Eq.[3]we have subtracted the T-independent 
susceptibility \w from xz(T) and xwB from M(H), cor- 
respondingly. For sample S2 (Fig. HI a scaling was ob- 
served for A G «0.5(1) with a scaling moment // ^3.5/i&. 
According to Ref. 20 the obtained moment can be related 
to a typical cluster size in crystal S2. The SH data are also 
in agreement with the G-like scenario (see below). There- 
fore, we ascribe the anomalous power-law of x{T) at low 
T < T* w 50 K to the formation of a quantum G-phase. 

2.2 Specific heat measurements 

Specific heat in the normal state We have found 
that the glassy magnetic subsystems (observed in the mag- 
netization measurements) do also contribute to the SH 
shown in Fig. [5^. In case of SG or CG phases the mag- 



netic contribution Ccg to the SH varies almost linearly 
at T < Tf like the usual electronic contribution in case 
of a FL [ 7 6 2 1 22 1 . Empirically, this behavior can be 
approximated by 



Ccg ~ IcgT + sqg2T 



or 



Ccg ~ ecGi.bT 



1.5 



(6) 



(7) 



where 7cc, £CG2 and ecgi.5 are CG constants. The ecgi.5 
contribution can be interpreted as originating from short- 
range 3D ferromagnetic (FM) spin waves which can exist 
in a FM clusters 11231 . also a linear contribution to SH can 
be expected for 2D FM spin waves. Then, the normal state 
SH of sample S 1 reads 



Cf (T) = 7cl T + C GG + (3 3 T 3 + p 5 T 5 



(8) 



where C GG is given by Eq.|6]and Eq.|7] 7^ is an intrinsic 
electronic contribution and ^3, f3§ are a lattice contribution. 
In case of a G-phase (sample S2), C{T) G /T oc x(T) is 
expected [24,25). Hence, for the SH we have: 



C P (T) S2 = lcl T + 7G T AG + /3 3 T 3 + f3 5 T ti 



(9) 



where A G pa 0.5(1) according to our magnetization data. 
To reduce the number of fitting parameters in Eqs. [8] and [9] 
we analyzed the difference: 



K^T) - C$\T)]/T = C* G - 1G T 



iAo-1 



(10) 



This allow us to exclude the lattice contributions /3 3 , (3$, 
as well as the linear electronic term j e \ which are all 
supposed to be nearly the same for both crystals, re- 
spectively. The fit of the experimental data by Eq. [10] 
is shown in Fig. [5p. i) In the case of Eq. [6] it gives: 
7 gi, = 36 mJ/mol-K 2 , e CG2 = 2.0 mJ/mol-K 3 and 
7c = 104 mJ/mol-K 1 ' 5 , respectively. Then, using in 
Eq. [9] obtained magnetic contribution we have estimated 
the intrinsic 7^ « 52 mJ/mol-K 2 for sample S2 with 
the lattice terms (3 % 3 = 0.59 mJ/mol-K 4 and j3\ = 
1.20- 10~ 3 mJ/mol-K 6 , respectively. The obtained /3 3 cor- 
responds to a Debye temperature Q % D 254 K. ii) In the 
case of validity of Eq.Qit gives: ecci. 5 = 14.9 mJ/mol-K 3 
and 7c = 75.3 mJ/mol-K 15 , respectively. Then the intrin- 
sic 7*1 f» 68 mJ/mol-K 2 is the same for SI and S2 crystals 
with slightly different lattice terms as compared to those 
obtained in the analysis (ii): /3" = 0.46 mJ/mol-K 4 and 
/?5* = 1.86- 10~ 3 mJ/mol-K 6 , respectively. This value of 
/3 3 l corresponds to a Debye temperature <9g ps276K. Both 
analysis give reasonable values of the lattice contribution 
(for example, in the case of Bao.6sKo.32Fe2As2 a Debye 
temperature of 277K was estimated l26l ) and essentially 
reduced 7e ] « 52 — 68 mJ/mol-K 2 as compared to nom- 
inal values 7 n 100 mJ/mol-K 2 l27ll8l obtained with 
out accounting for magnetic contributions. The SH data at 
B = 9 T shown in Fig. [5^ can be considered as a support 
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samples (see inset Fig.|5^). The fit below IK gives a resid- 
ual contribution for crystal SI 7£i-, «43 mJ/mol-K 2 and 
about 7^1 



to 7qq = 36 mJ/mol-K 2 estimated for the normal state 
using analysis (i). The closeness of these values would in- 
dicate that 7gQ is weakly effected by the SC transition and 
also excludes essentially a non-superconducting volume 
fraction for our the samples. The latter is also supported 
by the large absolute value of the SH jump at T c compared 
to the reported in the literature values f8l [27l . In the case of 
JSl it was observed that C p /T of the investigated crystals 
tends to zero at T = after AFM type magnetic transi- 
tion at T ~0.7K. This demonstrates that almost all itin- 
erant quasi-particals are gapped at T = 0. Therefore, we 
conclude that the large residual 7 res in the SC state of our 
samples is mainly due to the magnetic contribution from 
a CG. On the other hand, using ecgi.s = 14.9 mJ/mol-K 3 
from analysis (ii), we get 7^9 ~36 mJ/mol-K 2 . This value 



=46 mJ/mol-K 2 for S2 J29J. The 7^ is close 



Ires2 

68 mJ/mol-K 2 . 



In contrast to 



10 



Figure 5 (Color online) (a) The specific heat of our two 
K122 samples. Inset: low temperature part of the SH at 
zero field, (b) The plot [C^{T) - C^ 2 (T)}/T. 



is nearly a half of 7*] 
the conclusion obtained from analysis (i), it would mean 
that the CG phase in SC state is different from the CG in 
the normal state, since we exclude a large non-SC part of 
our samples. This can be possible, since itinerant electrons 
responsible for the RKKY interaction are affected by the 
SC transition. Thus, on this stage we cannot decide which 
analysis (i) or (ii) is more sophisticated. Therefore, we es- 
timate the intrinsic 7^ « 52 — 68 mJ/mol-K 2 . A more 
detailed report of the superconducting properties including 
microscopic considerations will be given elsewhere. 

2.3 Possible disorder induced quantum phase 
transitions Up to now the structural investigation of 
the cleaved surface of the samples such as EDX, XRD 
and SEM did not reveal any secondary phases ll9l fT0l . 



of analysis (ii), since this analysis provides essentially a 
better agreement between the experimental data and the 
fitting curves at low-T. However, we cannot exclude that 
large field enhances FM order in S 1 crystals and actually 
change the entropy of SG at low temperatures. 

Below Td ~6 K the data for S2 deviate from the fit- 
curve (Fig. [5]). At T C g ~4 K, slightly above T c , (see 
Fig. |5h) another magnetic anomaly is well visible in the 
SH data. Additionally, slightly above Tqg we observed a 
plateau-like feature in x{T) at relatively low fields (see 
Fig. 2k). We ascribe Td to the freezing of large cluster dy- 
namics in accord with the behavior expected for a quantum 
G-phase (see Ref. l28l ) followed by the final formation of 
a CG phase due to the RKKY interaction between the clus- 
ters at Tcg in crystal S2, too (for an illustration see also 
Fig.®. 

Specific heat in the superconducting state Mea- 
surements in SC state have shown that there is a large resid- 
ual Sommerfeld coefficient observed for all investigated 




concentration of MM (c„„) 

y MM 7 



Figure 6 (Color online) Schematic phase diagram of ex- 
trinsic magnetic moments (MM) driven quantum phase 
transitions. Notation of phases: SC - superconductivity, G 
- Griffiths, CG - cluster glass, SI, S2 - samples from this 
work. 



Copyright line will be provided by the publisher 



6 



V. Grinenko et al.: Disorded magnetism in superconducting KFe2As 2 



Therefore, we enforced to adopt a 'point' -defect model 
such as vacancies or interstitials of Fe atoms. To com- 
pare the amount of magnetic clusters contributing to a 
glassy phases of our samples, we calculated the magnetic 
entropy S m = §(C m /T)dT using the obtained above 
magnetic contributions. 5 T ^ 2 for crystal S2 related to CG 
and G phases between and T* m 50 K (where the quan- 
tum Griffiths behavior appears in the magnetization data 
Fig. is -0.074RJ/mol-Fe-K 2 . The estimate for crystal 
S2 related to the CG phases below Tf sw 87 K gives an 
essentially higher value than S^ 1 . In the case of validity of 
analysis (i) ^0.64RJ/mol-FeK 2 and for the case of analy- 
sis (ii) ~0.48RJ/mol-Fe-K 2 for the magnetic contribution 
have been obtained, respectively. Hence, we conclude that 
in the normal state crystal SI, it can contain up to 10 times 
more magnetic clusters than S2 does. Summarizing our ex- 
perimental observations of disordered magnetic phases in 
K122, we can propose a quantum phase transition of spin 
glass type with strong quantum G-phase effects (see Fig.|6]l 
driven by some tuning parameter p which is responsible 
for the formation of magnetic moments (MM) in K122. 
The physical nature of p should be unraveled in future in- 
vestigations such as spin and nuclear magnetic resonance 
and/or Mossbauer spectroscopy. These techniques can be 
helpful to estimate the amount and the distribution of MM 
in K122 single crystals. 

3 Conclusions To summarize, analyzing magnetiza- 
tion and specific heat data, we found out that even in high- 
quality KFe2As2 single crystals glassy magnetic behavior 
like in spin- , cluster-glasses and Griffiths phases may oc- 
cur near superconductivity and coexist with it. The mag- 
netic contribution is responsible for a large value of the 
nominal Sommerfeld coefficient 7„ ~100 mJ/mol-K 2 of 
this compound. The analysis of the SH data has shown 
that magnetic contribution amounts up to 50% of j n . In 
this way, the intrinsic value of the Sommerfeld coefficient 
7d ~ 52 — 68 mJ/mol-K 2 was estimated. We observed 
that various samples exhibit different disordered magnetic 
contributions depending on the amount and distribution of 
magnetic moments (MM). This suggests an extrinsic ori- 
gin of MM which can be caused by point defects such as 
vacancies or Fe interstitials. Therefore, we proposed a sce- 
nario of disorder induced spin glass type quantum phase 
transitions accomplished by strong quantum Griffiths ef- 
fects. Further investigations are required to elucidate the 
physical origin and the distribution of such MM. 
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